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Scaling Property of the F-AF Spin Chain Near the Exactly Solvable Point 
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We investigate the ground state of the J\-Ji spin-1/2 chain with J\ < and J2 > in 
the case that the nearest-neighbor J\ interaction in the z-direction has a weak anisotropy 
as Ji(l — a). We perform a perturbational analysis for small a and A = J 2 — |«7i|/4 with 
the exact solution of the unperturbed ground state for a — A = 0. The scaling property of 
the ground state energy is examined in detail. By the numerical diagonalization analysis of 
finite size systems, we found the phase boundary equation between the spin fluid and dirtier 
phases as a c = 14 Ac^ 3 . 
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1. Introduction 

Low-dimensional quantum spin systems have attracted great attention for many years. 
Among them, quantum spin chains with competing interactions have been fascinating sub- 
jects. This is mainly because these systems exhibit rich varieties of exotic ground states and 
phenomena owing to the quantum fluctuation enhanced by geometrical frustrations. For ex- 
ample, spontaneous dimerization, 1 ) quantum chiral phases, 2 ~ n ) l/3-plateau, 12 ~ 15 ) and singlet 
cluster solid 16 ' 17 ) have been reported. 

In theoretical studies of quantum spin systems, exact solutions are useful for constructing 
physical pictures of the systems. However, it is generally difficult to find an exact ground state 
of a frustrated quantum spin system. Despite this, exact ground states have been found in 
several models. For example, the exact dimer ground state for the antiferromagnetic J1-J2 
chain 1 ) is well known. Other examples are exact spin-cluster ground states for the pure and 
mixed diamond chains. 18 ~ 20 ) Even in the cases that the exact ground states are found, there 
still remains the difficulty to calculate correlation functions using them. 

Here, we are concerned with the F-AF chain depicted in Fig. 1 as another example that 
the exact ground state is known. This chain is a J\-Ji spin-1/2 chain with ferromagnetic (F) 
nearest-neighbor (NN) and antiferromagnetic (AF) next-nearest-neighbor (NNN) interactions 
(J\ < and J2 > 0). The ground state is fully ferromagnetic for J2 < | Ji|/4, 21 ~ 23 ) while it is 
a singlet state for J2 > |Ji|/4. At the point of J2 = |</i|/4, the explicit formulae of ground 



*Present address: Department of Physics, College of Humanities and Sciences, Nihon University, Setagaya-ku, 
Tokyo 156-8550, JAPAN 



1/15 



J. Phys. Soc. Jpn. 



Full Paper 




N-2 




JV-3 N-l 



Fig. 1. F-AF chain with the special open boundary condition. 




Fig. 2. Schematic phase diagram of the F-AF chain. The exactly solvable point is indicated by the 
dot at a = A = 0, and the phase boundary between the spin fluid and dimer phases is shown by 
the bold solid line. 



states have been found and studied. 24 ~ 27 ) In particular, Hamada et al. first found an exact 
solution in a resonating valence bond form. 24 - 1 Recently, we have reported all the degenerate 
ground states in the explicit formulae and further proposed the recursion formulae of the exact 
solutions. 27 ^ These exact solutions are extended to include the case with bond alternations in 
the NN and NNN interactions. The greatest merit of the recursion formulae is that physical 
quantities such as correlation functions can be calculated using them. 

The ground state near the exactly solvable point has also been examined. Several authors 
studied the issue by adding two deviation terms to the Hamiltonian at the exactly solvable 
point. 28 ~ 31 ) One is the anisotropy term with the anisotropy parameter a, which displaces 
the exchange energy of the z-direction as Ji(l — a). The other is the NNN interaction term 
measured from the exactly solvable point; namely, the NNN coupling constant is the parameter 
A = J2 — I Ji |/4. The numerical studies show that the ground state for a > and A = belongs 
to the spin fluid (Tomonaga-Luttinger liquid) phase and that for a = and A > belongs 
to the dimer phase. 28_3 °) For a, A <C 1, Dmitriev and Krivnov 31 ) discussed a scaling property 

of the ground state energy. They also argued about the phase boundary between the spin 

4/3 

fluid and dimer phases and found the equation of the phase boundary as a c ~ 13.9A C using 
the scaling property of the ground state energy in the dimer phase. The sketch of the phase 
boundary is shown in Fig. 2. 

In this paper, we investigate the ground state behavior of the anisotropic F-AF chain near 
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the exactly solvable point. We consider the Hamiltonian of the F-AF chain at the exactly 
solvable point of a = A = as an unperturbed system and the other terms mentioned above 
as perturbations. Using the recursion formulae of the exact solutions of the ground states, 
the ground state energy can be calculated within the first order term of the perturbation. 
From the result of the first order perturbation, we derive a scaling form of the ground state 
energy in the spin fluid phase. This scaling form is confirmed numerically using the exact 
diagonalization. The detailed expression of the phase boundary between the spin fluid and 
dimer phases is determined from the scaling property with numerical data. 

This paper is organized as follows. In §2, the Hamiltonian of the unperturbed system and 
the perturbation terms are explained. In §3, the exact solution of the ground state of the 
unperturbed system is given. In §4, within the first order perturbation theory, we discuss the 
phase boundary between the spin fluid and dimer phases. In §5, the detailed expression of the 
phase boundary is derived using the scaling property of the ground state energy. Finally, §6 
devoted to summary and discussions. 

2. F-AF Chain with Perturbations 

We examine the effects of perturbations near the exact solutions for the F-AF chain at 
the exactly solvable point. The total Hamiltonian consists of the unperturbed Hamiltonian 
Tio and the perturbation terms V\ and V a : 

n = n + v x + v a , (i) 

"Ho = ^2 (-SnSn+l + jS n S n+2 J , (2) 
V X = \J2(*nSn+2- J), (3) 

V a =aJ2 s nSn+l, (4) 
n 

where s n is the spin- 1/2 operator at the n-th site, and the summations are taken over the 
total number N of spin sites. We used the energy unit of J\ = — 1, so that eqs. (2) and (3) 
indicate J 2 = A + 1/4. 

To treat the exact solutions for the unperturbed Hamiltonian Hq at the exactly solvable 
point, we write it as a sum of local Hamiltonians: %q = J2 n h n , with 

K = --^(SnSn+l + S n+1 S n+2 ) + -S n S n+2 . (5) 

Here, the system size N is assumed to be an even number. Owing to the form of the Hamil- 
tonian, it is possible to use a special open boundary condition (OBC) with n = 1 to N — 2, 
shown in Fig. 1, or a periodic boundary condition (PBC) with n = 1 to N. The ground states 
under the OBC are degenerate with respect to all the magnitudes and z-components of the 
total spin. We have found the exact recursion formulae for all the degenerate ground states. 2? ) 
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We choose the OBC, since we need all the degenerate ground states to perform first-order 
perturbation calculations. The degeneracy of these ground states of Hq is resolved by the 
perturbation terms V a and V\. 

The perturbation term V\ changes the NNN interaction parameter J2 from 1/4 at the 
exactly solvable point to A + 1/4. Accordingly, for A < 0, the system shifts from the exactly 
solvable point into the ferromagnetic phase, if a = 0. In contrast, for A > 0, it shifts into the 
dimer phase. 

The other perturbation term V a represents an anisotropy of the NN interaction. If V a is 
Ising-like (a < 0) and A = 0, the system shifts into the ferromagnetic phase where the total 
spin direction is along the z-axis. In contrast, if it is XY-like (a > 0), the system shifts into 
the spin fluid phase. In this case, the ground state is approximated by a linear combination 
of the degenerate ground states of T~Lq in the first-order perturbation. We concentrate on this 
interesting case of a > 0. Also, we consider the anisotropy only for the NN interaction and 
not for the NNN interaction for simplicity. The latter brings no drastic change as will be 
discussed. 

3. Exact Solution for a = X = 

We summarize the exact recursion formula for the unperturbed Hamiltonian T~Lq with the 
OBC, which describes the uniform and isotropic F-AF chain at the exactly solvable point. 
The general formula in the case with full bond alternations in the NN and NNN interactions 
is given in our precedent paper. 27 ) 

Let \j, m)jy be the ground state of the TV-site chain for S to t = j and S% ot = m, where 
S'tot and S* ot are the magnitude and the z-component of the total spin Stot = Yln=i s ™' 
respectively. The ground state \ j, m) 7V+2 of the (N + 2)-site chain is formed by incorporating 
spin states of two extra 1/2 spins at the (N + l)-th and (N + 2)-th sites. The triplet and 
singlet states for the extra spins are expressed as 

|tl) = I tiV+ltiV+2); |t-l) = I lN+llN+2), 

I *o) = ^=(| tN+llN+2} + I J-JV+ltjV+2)), 

I s ) = ^=(1 tiV+l4iV+2) - I iN+ltN+2))- (6) 
Then \ j, m) at +2 is written as 

\j, m) N+2 =a N (j)^2,C(\,n)\j + l,m- fi) N ® |t M ) 



+bN(j)^2C(0,fi)\j,m - n) N ® \tfjt) 
+CN{j)^2C(-l,n)\j -l,m- n) N (8) |t M ) 
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+dN{j)\j, m - n) N ® |s) (7) 

with the Clebsch-Gordan coefficients 

C(i/,/i) = j + is,m- n\j,m) (8) 

for fi = — 1, 0, and 1. The coefficients {a^, bjy, cn, d^} are independent of m because of 
the rotational symmetry of T~Lq. If an unreasonable state like \(N/2) + l,m)jv appears in the 
right-hand side of eq. (7), we regard the coefficient for it as zero. 

We impose the condition that eq. (7) is a common ground state of the local Hamiltonians 
/ijv-i and h jy . Then we obtain the recursion relations for {a^, bw, cn,cIn} with respect to the 
system size N as 

clnU) _ a N _ 2 (j) 



Mi) 



djv- 2 (i + i)' 

bN-2{j) 



dN{j) d N - 2 {jY 
cjv(i) C N - 2 (j) 



dN(j) d N - 2 (j-l)' 
Also, the following relations for the same iV are satisfied: 



dwti) 



lj - 1 b N (j - 1) 



j + 1 d N (j - 1) 



+ 44 



i + i 



(9) 
(10) 

(11) 

(12) 



and 



dN{j) 



y/j^lb N (j - 1) + ^d N (j - 1 
y/j=Tb N (j - 1) + 3^jd N (j - 1) 



(13) 



(14) 



V(2j + l)(2j-l)ajv(j - l)d N (j - 1) 
From the normalization condition of \j, m)jv+2j w e have 

«jv(j) 2 + Mi) 2 + cjv(i) 2 + d N {jf = i. 

The recursion relations eqs. (11)-(14) determine the coefficients {a^, b]y, cn, d^} by start- 
ing from the initial values for N = 4 except for coefficients with special values of j. Since 
the denominators of the right-hand sides of eq. (9) for j = (N/2) — 1 and of eq. (10) for 
j = N/2 become zero, we separately evaluate b]\r(N/2) and aAr((A r /2) — 1) using eqs. (12) and 
(13), respectively. The initial values of the coefficients and initial ground states are given in 
Appendix. 

Using the above recursion relation with respect to the system size N, the matrix element 
(j' ,m'\Q\j,m) n for an arbitrary Q can be decomposed into the summation of the two-spin 
matrix element. Thus, we can easily calculate an arbitrary correlation function for large N. 



5/15 



J. Phys. Soc. Jpn. 



Full Paper 



0.88 - 




0.0 0.5 1.0 1.5 

a/X 

Fig. 3. Ground energy and first excitation energy for N = 1 8 and 20. 

4. First Order Perturbation 

The perturbation term V a + V\ in eq. (1) resolves the degeneracy of the ground states 
of Hq with the OBC. We have considered the case of even N for simplicity. We assume that 
the ground state of % exists in the subspace of S^ ot = 0. Then, the zeroth order ground state 
\ip0)N is written as a linear combination of the basis vectors of the subspace: 

N/2 

\*Po)N = ^2b j \j,0) N . (15) 

j=0 

The coefficient vector b = (bo, b\, ■ ■ ■ , frjvya) is the eigenvector for the smallest eigenvalue of 
the matrix v whose (j,f) element is defined as Vjj> = (j,0\V a + V\\j',0)N- The first excited 
state \iPi)n is also written in the same form as eq. (15) but with the coefficient vector for the 
next smallest eigenvalue. The shift of the ground state energy for the first order perturbation 
is given by 

8E$(N) = (ipo\(V a + V x )\^o)n, (16) 

which is just the eigenvalue of v. 

The matrix element Vjji is nonzero only for j = j' or j = j' ± 2. Hence, we have b2 n = 
for all n, otherwise &2n+i = for all re. Therefore, at least either bo or 6jv/2 is zero if N/2 is 
odd, while both of them can be nonzero if N/2 is even. For a = and A > (dimer phase), 
eq. (15) becomes \iPq)n = |0, 0}jv- Hence, the state \ipo)N for the dimer phase contains the 
state |0,0)jV) namely bo ^ 0. On the other hand, for a > and A = (spin fluid phase), 
eq. (15) becomes \ipo)N = \N/2,0)n- Hence, the state |^o)iV f° r the spin fluid phase contains 
the ferromagnetic state \N/2, Q)n, namely 6jy/2 7^ 0- Thus, for odd N/2, the eigenvector for 
the dimer phase and that for the spin fluid phase are orthogonal, so that the energy level 
crosses at the phase boundary. We illustrate the level crossing for N = 18 in Fig. 3. This fact 
is observed not only in the first order perturbation but also in the numerical diagonalization 
of the total Hamiltonian with the OBC. This is an effect of the OBC for a finite size system. 
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The crossing points a c for N < 90 are plotted in Fig. 4. The solid line is obtained by 
fitting using data for N = 62 ~ 90. We obtain a c — > in the thermodynamic limit N —> oo. 
Therefore, a c can be expressed as a function of A in which the exponent of A is larger than 1. 
This is consistent with the phase boundary given in ref. 31 and the following section. 

5. Scaling Property and Phase Boundary 

We argue about the scaling property of the shift of the ground state energy for small a 
and A from the ground state energy Eq of the unperturbed system %q. Before examining the 
general case correctly, we first treat the special case of A = under the PBC, following the 
argument of Dmitriev and Krivnov. 31 ) Based on the perturbation argument, they introduced 
a scaling parameter 

(17) 

h, k — Vjq 

where Ek is a typical low eigenenergy of T~Lq and \ip k ) is the excited state belonging to E k . 
In the most divergent contribution, the ground state energy and the corresponding eigenstate 
are respectively expanded in powers of x as 

5E (a) = E (a) - E = (Vol Wo) Yl c ^> ( 18 ) 

i 

\M<*)) =J2dkX 9k Wk), (19) 
k 

where Cj's and d^s are coefficients and q^s are non-negative integers. The 1-magnon excitation 
spectrum of Ho is written as 32 ) 

E k -E = (l- cosfc) - ~(1 - cos2/c) ~ ^ ~ -ij. (20) 

The numerator (tpj\V a \ipk) of x is analyzed and found to be independent of N for large iV: 31 ) 

(V>;IWfc)~a. (21) 
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Thus, the scaling parameter x of the shift of the ground state energy is given by 

x = aN A . (22) 

Since the energy is an extensive variable, eq. (18) can be written using the scaling function 
f(x) as 

5E Q (a) = Na b/4 f{x). (23) 

This equation is confirmed by the numerical diagonalization with the PBC. 31 ^ 

To use the scaling form eq. (23) in our calculations, we need finite size correction in the 
case of the OBC. Using the exact recursion relations presented in the last section, we calculate 
5Eq{q) in the first order perturbation with respect to a. The result is shown in Fig. 5. The 
energy shift is represented clS cl function of N as 

6eW{N) = lim = »k 

a-to a a 

The iV-dependent term in parentheses is the correction from the OBC, which is not negligible 
for iV < 24. Therefore, the corrected energy shift has the following scaling form 

^9 = Na^foMx). (25) 



0.509 x N 0A95 

The scaling function obtained from the numerical diagonalization result for N = 16 ~ 24 with 
the OBC is shown in Fig. 6. 

Now, we investigate the general case with both the perturbation terms, V a and V\. Here, 
we assume that the parameter A takes a sufficiently small value so that the ground state 
belongs to the spin fluid phase. The scaling parameter y caused by V\ is introduced as 

v ^mm. (26 ) 

Hih — -C/n 
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The iV-dependence of the matrix element {ipj\V\\i/jk) is written as 

(tf'ilWk) ~ AiV e (27) 
with an exponent 0, which will be determined later. Then the scaling parameter is represented 



as 



XN e+i . 



y = AiV" 1 \ (28) 

We do not adopt the formula y = XN 2 that is used in ref. 31, since we find no adequate 
argument to justify it. Using the scaling function f(x,y), the ground state energy is written 
as 



5E (a, A) = E (a, A) - E = Na 5 / 4 f(x, y). 



(29) 



To determine the exponent 6, we consider the energy correction in the first order of A. 
Using eqs. (19) and (27), we have 

E (a,X)-E (a,0) = {M<x)\V*\M<*)) 

jk 

= XN e g(x), (30) 

where g(x) is a scaling function depending only on x. We cannot evaluate (ipo(a)\Vx\ipo(a)} 
directly in the first order perturbation with respect to a, since we have no exact solutions 
for excited states to calculate the O(a) term of \ipo(a)}. Therefore, we separately calculate 
Eq(oc, A) and Eq(oc, 0) in the first order of a; the corrections are written as 

E (a, 0)=E + a SE® (JV) + 0(a 2 ), 



E (a, X)=E + a 6E$ a (N) + 0(a 2 ). 



(31) 
(32) 
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Fig. 7. A dependence of the ground state energy in the first order perturbation. 
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Fig. 8. Size dependence of the coefficient e(N). 



The former correction 5Ea\N) has been given in eq. (24). We estimate the latter correction 



A 

a 



s n s n+l + I s n s n+2 . 



(33) 



using the recursion relations in the last section. The differences of the corrections for several 
values of N are shown in Fig. 7. The result is represented in the form 



5E^(N)-5E^(N) = e(N)^. 



(34) 



The function e(N) is linear for large N as is shown in Fig. 8. Consequently, we have 

(Vo(a)IWo(a)) ~ XN-\ (35) 



and then the scaling parameter y is given by 

AiV^ 1 



y 



(36) 



10/15 



J. Phys. Soc. Jpn. 



Full Paper 



0.(10 




A 
□ 


N=14 
N=16 
JV=18 


"U^ 

\ 




• 


N=20 








N=22 


o 




■ 


N=24 *x 


a 




1 /Na 


1/4 = 0.2 




0. 



1.0 1.5 2.0 2.5 3.0 

Fig. 9. Scaling function f(x, y) vs y — NX 1 / 3 for a fixed value of 0.2 for x = 1/ (Na 1 / 4 ). The crossing 
point (x c , y c ) corresponds to the phase boundary for x c = 0.2. 



We confirm the scaling form eq. (29) by the numerical diagonalization under the PBC for 
finite size systems with N = 14 ~ 24. The y 1 / 3 dependence of the scaling function f(x,y) 
with a fixed x for the ground state energy is shown in Fig. 9. Clearly, all the data lie well on 
a unique curve. We also show the same for the first excitation energy 5Ei(a, A). All the data 
also lie well on a curve. This means that the low excitation energy is also described by the 
same scaling parameters x and y. 

Extending the present scaling form for larger A, we find a well-defined level crossing point 
y c for each x c . This point corresponds to the phase boundary between the spin fluid and dimer 
phases. At this point, eliminating N from x c = a c N 4 and y c = A c iV 3 , we have the following 
equation, 

«c = -%\ A J 3 . (37) 

Here, y c is estimated by the numerical diagonalization for each fixed x c . The x c dependence 
of x c /yt^ 3 for N = 22 is shown in Fig. 10. In the thermodynamic limit (x = aN 4 — > oo), we 
have x c /yt^ — V 14 ± 0.1 and finally arrive at the following expression of the phase boundary: 

a c ~ 14A^ 3 . (38) 

6. Summary and Discussion 

To summarize, we studied the phase boundary of the F-AF spin chain with an anisotropy 
and an extra NNN interaction, which are characterized by the parameters a and A = J2 — 
I Ji|/4, respectively. The analysis is perturbational with the exact solutions in the unperturbed 
case of a = A = 0. We investigated the scaling property of the ground state energy in the spin 
fluid phase. We numerically confirmed the scaling form and arrived at the phase boundary 
equation between the spin fluid and dimer phases. 

Here, we mention the following two points in comparison with our results and the result 
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Fig. 10. x c dependence of x c /yt for N = 22. 

of a previous study by Dmitriev and Krivnov. 31 ) First, the scaling parameter y obtained 
by Dmitriev and Krivnov (y = XN 2 ) has a different expression from eq. (36), because they 
adopted only the ferromagnetic ground state as a unperturbed state. The scaling parameter 
of the ground state energy in the spin fluid state is correctly described by eq. (36). Second, 
our result of the phase boundary has the same form as the result they obtained. We used a 
derivation based on the perturbation in the spin fluid phase, while they used the perturbation 
in the dimer phase. Thus, eq. (38) is the definitive formula of the phase boundary between 
the spin fluid and dimer phases. 

When the NNN interaction has the same anisotropy as that of the NN interaction with an 
XY-type, the similar scaling property discussed in §5 can be obtained. However, the scaling 
function plotted using the numerical diagonalzation for N = 16 ~ 20 does not completely lie 
on one curve. If we use larger systems, the diagonalization data will lie on one curve. The 
phase boundary estimated by the chain for N = 22 is a c ~ 15Ac^ 3 . 

In the case of the Ising-type anisotropy of the NN interaction (a < 0), the first order 
phase transition between the dimer and ferromagnetic phases exists. However, the detailed 
investigation using a finite size calculation is disturbed by the complicated size dependence 33 ) 
of the ground state energy in the dimer phase. The meanfield approximation predicts that 
the phase transition takes place at a c ~ — A^ 3 . 31 -* Within the first order perturbation used in 
§4, we obtain the result that the linear term of a c with respect to A c becomes zero, which is 
consistent with the meanfield result. In the case that the NNN interaction also has the same 
anisotropy as that of the NN interaction with an Ising-type, there appears an intermediate 
state between the dimer state and the fully polarized ferromagnetic state. 34 ) Unfortunately, 
this intermediate phase does not appear by the first order perturbation. The boundary between 
these phases can be roughly estimated by the numerical diagonalization, however, the complex 
size dependence disturbs the extrapolation of the system size N. Detailed study in this regime 
remains in the future. 
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Tabic A-l. Initial coefficients {04,64,04,^4} for the recursion relations eqs. (9)-(ll). 
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Appendix: Initial state of the recursion relation 

The smallest chain decomposed to the local Hamiltonian eq. (5) is N = 4 for even N. 
Thus, the initial ground states for eq. (7) are given as follows: 

I tttt), 
(1,1, 1,1) -«i 



2,2)4 = 
2,1)4 = 
1,1)4 = 
2,0)4 = 
1,0)4 = 

0. 0)4 = 

2,-1)4 

1, -1)4 



(3,1,-1,-3) -iti 
2^5 

(1,1,1,1,1,1) - up 
v/6 

(2,1,0,0,-1,-2) - wo 



10 

(3,2,-3,-3,2,3) - up 

(UMHm 
2 

_ (3,1,-1,-3) ju-i 
2^5 



(A-l) 



\2, -2) 4 = 1 1114), 

where u, = (| tttt), I tttt), I tttt), I tttt))*, «o = (| tttt), I tttt), I tttt)), I tttt), I tttt 
), I tttt))* and u 1 = (| tttt), I tttt), I tttt), I tttt))*- The initial coefficients {04,64,04,^4} 
for the recursion relations eqs. (9)-(ll) are given in Table A-l. Thus, the ground state \ j,rn)N 
with arbitrary even N is constructed recursively by eqs. (7) and (9)-(14). All the coefficients 
can be determined as positive values by a suitable choice of the sign of each ground state. 
Thus, the coefficients {a^, b^, cat, djy} are determined uniquely. Therefore, the ground state 
in the sector of fixed S to t = j and S* ot = m is nondegenerate. The total degeneracy of the 
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ground states is (N + 2) 2 /4 for even N. 
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